Fluctuations in finite critical and turbulent systems 
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We show that hyperscaling and finite-size scaling imply that the probability distribution of the 
order parameter in finite size critical systems exhibit data collapse. We consider the examples of 
equilibrium critical systems, and a statistical model of ecology. We explain recent observations that 
the probability distribution of turbulent power fluctuations in closed flows is the same as that of the 
harmonic 2DXY model. 
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Critical systems with identical symmetries, dimension 
and exponents are denned to be members of the same 
universality class; but must they also share the same 
probability distribution for the fluctuating variables or 
order parameters? And if two systems do indeed exhibit 
the same probability distribution, are they necessarily in 
the same universality class, in the conventional meaning 
of the term? 

Recently, light has been shed on these issues by studies 
of the probability distribution functions (PDFs) of fluc- 
tuating quantities in finite critical systems. Bramwell, 
Holdsworth and Pinton (BHP) [Q observed data col- 
lapse for the two-dimensional XY model (2DXY) in the 
spin-wave regime at low temperatures and in statistical 
models of nonequilibrium dynamics which exhibit self- 
organized criticality (SOC). Even more remarkably, they 
found that the power fluctuations in a closed turbulent 
flow exhibit exactly the same form of data collapse, with 
a scaling function that is indistinguishable from the afore- 
mentioned statistical critical models. Taken at face value, 
these observations lend support to the notion that fi- 
nite Reynolds number (Re) turbulence is indeed a critical 
state, and that there is a kind of super-universality be- 
tween systems with different dynamics and even dimen- 
sionality. 

The purpose of this Letter is to show that finite-size 
systems that are in the critical regime should be expected 
to exhibit just this sort of data collapse. The system in 
question may be either an equilibrium system near its 
critical point, or a nonequilibrium system that attains a 
critical state through fine-tuning or other mechanism for 
achieving scale invariance. However, we are unaware of 
any reason a priori to expect that the probability distri- 
bution should be super-universal, and indeed we exhibit 
a counter example. Finally, we argue that the apparent 
agreement between magnetic systems and experiments on 
closed turbulent flows, while interesting and genuine, is 
not indicative of the intrinsic behavior of turbulence; we 
propose an explanation for the observations that appears 
to explain not only the data collapse but the Reynolds 
number dependence as well. 



Probability Distribution Data Collapse:- Let us now re- 
view in more detail the findings of BHP. They examined 
the 2DXY model in the harmonic approximation for tem- 
peratures well below the Kosterlitz-Thouless transition 
Tkt- Although in an infinite system the magnetization 
M should be identically zero, for a finite system there are 
large fluctuations, and they measured the probability dis- 
tribution function (PDF) for a range of system sizes and 
temperatures. 

The scaled PDFs of the magnetization below Tkt col- 
lapse onto each other for different system sizes and tem- 
peratures, provided one works in the harmonic approxi- 
mation and the correlation length, £, is larger than the 
system size. The scaling necessary to achieve this data 
collapse is that the independent variable (rj = \M\) is 
replaced by y = rj — (rj) /a, where a is a measure of the 
width of the PDF, such as the width at half maximum. 
A similar data collapse was seen in a statistical model 
of ecology [0, where the scaled probability distribution 
of species abundance in a region was found to be inde- 
pendent of its area. Pinton et al. Hflj performed experi- 
ments on confined turbulent flows maintained at constant 
Reynolds number and looked at the PDF of power fluctu- 
ations. These too showed data collapse across different 
Reynolds numbers. Remarkably the PDFs for the tur- 
bulence data and the 2DXY model overlap within the 
precision of the data. 

Subsequently a number of SOC systems, such as the 
Bak-Tang-Weisenfeld sand pile model ||, the Sneppen 
depinning model ||, the auto igniting forest fire model 
and a model for granular media [|| , have been studied 
and they too seem to show data collapse with a PDF very 
similar to that observed in the turbulence experiment 
||. The same holds true of the 2D Ising and 2D site 
percolation models as well. This seems to suggest that 
the scaling form is independent of system attributes such 
as symmetry (discrete or otherwise), state (equilibrium 
or otherwise) etc. 

A similar phenomena was noted a long time ago by 
Nicolaides and Bruce Jl0| , who were interested in the 
question of whether a universality class is defined by the 
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values of the critical exponents, or whether the probabil- 
ity distributions were common to members of the same 
universality class. They found that the PDF of the two 
dimensional Ising, spin 1, and <p A models all had the same 
form in finite systems. 

In fact, the issue of data collapse is linked with the 
existence of hyperscaling. Let us begin with a discussion 
of data collapse in equilibrium systems. In finite size 
magnetic systems, the PDF has a scaling form in the 
critical region 0, P L {M) = L p / u f{L/^ MlPl v ), where 

/ is a scaling function, /3, v are the critical exponents 
for the order parameter M and the correlation length £ 
respectively, and L denotes the linear dimension of the 
system, assumed to be in d dimensions. 

Near the critical point the correlation length becomes 
larger than the system size and making the unjustified 
assumption (a priori) that / is analytic in the first argu- 
ment, we obtain P L (M) = L^I V ]{ML^I V ). The sum rule 
for the static susceptibility \t implies that ksTxT ~ 
L a [(M 2 ) - (M) 2 }. Notice that the quantity in the brack- 
ets is the variance of the probability distribution of the 
magnetization. The finite size scaling form for the sus- 
ceptibility is xt = Lpl v F(L/£) In the limit xjL goes to 
infinity the scaling function F tends towards a constant. 
Thus the measure of the width of the fluctuations, a, 
scales as a 2 = (M 2 ) - (M) 2 ~ L^ v ~ d 

BHP found that aP(M/a) is a universal function, in- 
dependent of system size. To test this, compute the func- 
tion aP(M) = lpl 2v - d L^l l - , f(ML l v ), 

aP(M/a) = L^ +2 f 3 - d ^! 2v f{M/aLf } ! v+ ~<l' 2v - d / 2 ) (1) 

Combining the hyperscaling relation, 2 — a — dv and 
Rushbrooke scaling, a + 2(3 + 7 = 2 gives 7 + 2/3 = dv. 
Given this identity it follows that all L dependence in 
(1) disappears and we are left with a statement of data 
collapse: 



aP{M/a) = f(M/a) 



(2) 



Thus as long as finite-size scaling holds true near the 
critical point and hyperscaling is obeyed, the data, for 
different sizes, fall on top of each other for a given system. 
To our knowledge, this was first observed empirically by 
Nicolaides and Bruce [ fToj ] . 

There are interesting consequences for the moments of 
the PDF, P, illustrated here for the first two. The mean 
of the distribution is 



(M) = 



f MP{M)dM f zF(z)dz 



J P(M)dM J F(z)dz 



(3) 



where F is a scaling function 

The integral is a pure number, so that the ratio of the 
mean and the variance is independent of L ({M) ~ a). 

The moment relation is a direct result of hyperscaling, 

(M) ~Lf ) l v 1 (j~Lr<l 2v - d l 2 



(M) 



L 



-2f3+j~dis/2u 



(4) 



Hyperscaling then requires the L dependence in the last 
line above to vanish. This relation between the mean 
and variance has been explicitly calculated for the 2D 
XY in the spin wave approximation ||. The moments 
of the order parameter satisfy, fi n = 9n{g2/2)~ n / 2 cr n and 
~ Mi ~ c™ implying 



Mi ~ (M) 



(5) 



Probability Distribution Data Collapse in an Ecology 
Model:- This scaling is fundamental to data collapse and 
holds even for systems not in thermal equilibrium. In 
particular it holds for the ecology model mentioned ear- 
lier. Harte et al. £g| proposed a model for the species 
abundance distribution observed in nature. It has been 
empirically observed that the number of species S in a 
patch of area A obeys a scaling law, 5 ~ A z . Presumably 
this is an asymptotic in time form of a more general dy- 
namical model, as the time evolution within this model is 
not specified. Nevertheless one can construct a recursion 
recursion relation for the PDF, Pi(n), the probability for 
finding a species with n individuals in a patch i. Patches 
i are constructed by successive bifurcation of an initial 
biome. The key ingredient of the model is the assump- 
tion of self similarity, which forces the species-area law, 
but allows the probability distribution for the number of 
species resident in a biome to be calculated. If a species 
is found in area Ai, then there is a non-zero probability a, 
for finding it in one of the two halves (area Aj+i = Ai/2,) 
which is independent of i. The independence on scale 
i of the probability a gives rise (or more accurately, is 
equivalent) to the species-area rule, with a — 2~ z . The 
resulting equation turns out to be 

n—1 

P i {n)=xP i+ x{n) + {l-x) y £ j P i+1 {n-k)Pi + i{k) (6) 

fc=i 

where x = 2(1 — a). The variance of this PDF was com- 
puted by Banavar et al. |12|: 



m — i — l 

E 

j=0 



{2-x) j (2-x) 



m — i — l 



(7) 



where m is defined as the maximum number of times 
the system can be subdivided before no species exist on 
the smallest patch (P m (l) = 1). When the system size is 
much bigger than the size of this smallest patch, the vari- 
ance obeys a simple recursion relation, ai/a^i ~ 2 — x. 
The data collapse observed for this model is the state- 
ment o-jPi{n) = (jj+iPj+i(n') and njai = n 1 /oi+\. It can 
be shown [O that the PDF satisfies a finite-size scaling 
relation of the form Pj(n) = l/nf(n/Nf), where cf> is a 
crossover scaling exponent, Ni is the number of individ- 
uals of all species in biome patch i, and Nf +l = 2~^Nf; 
then the data collapse is equivalent to a relation between 
the exponents. 



/( 



(2 - x)Nf 2-* 



(8) 
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leading to 2 — x — 2^ which is equivalent to <p + z = 1. 

This relation is nothing but the hyperscaling rela- 
tion observed in the magnetic system. This can be 
seen by considering the asymptotic forms of the mean 
and variance of the PDF, < n > ~ and a ~ 

(2 - x) m ~ 2 m< ^ ~ N*. This immediately gives us the 
moment scaling which was crucial to data collapse and 
hyperscaling, a ~ (n). It is remarkable that the ecology 
model shows this behaviour for its moments - an un- 
forseen consequence of the assumption of self-similarity. 

Power fluctuations in a turbulent flow:- Having estab- 
lished the connection between hyperscaling (moment re- 
lationships) and data collapse let us turn our attention 
to the case of confined turbulent flows. The experiment 
consists of a closed cylinder in which a turbulent fluid is 
driven at the top and bottom by counter-rotating plates 
with vanes, moving at the same mean angular frequency. 
The PDFs for power fluctuations II were measured [^|^] 
for different Reynolds numbers (Re) and found to exhibit 
data collapse. However the ratio of the mean to the vari- 
ance depended on Re: n rms /n ~ Re~ a , a ~ 0.33. 

This observation shows that the Reynolds number, and 
hence the extent of the inertial range, is not the parame- 
ter that controls the nature of the fluctuations: it is not 
analogous to the finite size of the system, as has been 
suggested previously. If this suggestion were correct, hy- 
perscaling would have required that the ratio of mean to 
variance be independent of Reynolds number. 

Our analysis suggests that the relevant parameter, in 
the case of the confined flows, that controls the nature 
of the fluctuations is not the Reynolds number but the 
physical system size itself. In addition to the scaling 
analysis mentioned above, the reasons for this is two fold. 
First, the experiments were performed by changing either 
the fluid (i.e. viscosity) or the angular velocity (see Ref. 
3 and 4 for details of the experiment) and not by changing 
the spatial dimensions of the flow. Second, it was also 
noted that an open flow did not exhibit any non-Gaussian 
fluctuations in the total power dissipated. In other words 
the fluctuations are not governed by the degrees of the 
turbulent flow (i.e. Reynolds number) but the degrees of 
freedom of coherent structures that form in these system 
which in turn depends on the system size, L. 

We hypothesize that the flow is composed of a tur- 
bulent region around the top plate with a mean angular 
momentum, another oppositely directed turbulent region 
around the bottom plate, and an interfacial region be- 
tween them - a shear pancake. Experiments were per- 
formed in both open and closed geometries, which dif- 
fer by the absence or presence of a confining cylindrical 
wall. In the open geometry, the power fluctuations of 
each plate were non-Gaussian, and negatively skewed, 
while the total power (sum of the power measured at 
each plate) fluctuations were Gaussian. In the closed ge- 
ometry all of these three quantities were non-Gaussian, 
with a scaled probability distribution apparently indis- 
tinguishable from that of the 2D XY model simulations. 

In the open geometry, the shear pancake can shift its 



mean vertical position instantaneously as shear energy is 
dissipated horizontally, but in the closed geometry, this is 
not possible because of the walls. In the open geometry, 
as the shear pancake moves, the net turbulent energy 
in the upper half of the cell will increase/decrease while 
that in the lower half of the cell will decrease/increase: 
hence, we anticipate that the power fluctuations should 
be anti-correlated, in agreement with observations 

In the closed geometry, we can describe the shear pan- 
cake by its height h(x, y) above the x, y plane positioned 
parallel to and equidistant from the rotating plates. The 
shear pancake experiences random fluctuations from the 
turbulent flows, which drive it with an effective, Reynolds 
number dependent temperature T(Re) and Boltzmann 
distribution P h (h) cx exp[-S(h)/T]. The action S(h) 
describes local height fluctuations of the shear pancake, 
and the power dissipation of the system should be ex- 
pected to depend in some way on the friction between 
the two counter-rotating turbulent flows i.e. on the sur- 
face area of the shear pancake, and not its absolute mean 
height. Any local fluctuation increases the total surface 
area of the shear pancake and hence increases the power 
dissipation in both flows. Hence, in this case, we antici- 
pate that the power fluctuations would be correlated, in 
agreement with observations [^5) . Thus the action should 
be of the form 

S(h)(X J dxdy^/l + (Vh) 2 « / dxdy(Vh) 2 (9) 

which is the Hamiltonian of the 2D XY model in the spin 
wave approximation, making the identification of h with 
the phase 9. 

In order to complete the dictionary between the PDFs 
of the 2D XY model magnetization and the turbu- 
lent power, we note that the power dissipated is by 
hypothesis proportional to S{h} and thus is linearly 
dependent on (V/i) 2 . Similarly, the absolute value 
of the magnetization per spin M, measured in the 
2DXY simulations (for N spins Si) is given by M = 

^(Ef =1 £)-(Ef =1 ^-)/iv « i - ^E«(o< - ^) 2 > 

showing that the probability distribution of magnetiza- 
tion fluctuations is indeed the counterpart of the power 
fluctuations in the turbulent flow. It follows from our 
assumptions above that the probability distributions 
should be identical. Note that we are definitely predict- 
ing that the power probability distribution should be that 
of the 2DXY model in the harmonic approximation only. 
Indeed, we have simulated the 2DXY model for temper- 
atures near the Kosterlitz-Thouless transition, and find 
significant deviations from the probability distribution 
for the turbulent power and the 2DXY model in the har- 
monic approximation as shown in figure 1 (and has also 
been noticed by Palme et al. [fl6||). 

Now, we discuss the Reynolds number dependence of 
the ratio ixja. By hypothesis, this should be given by 
using the temperature dependence of \xja from the har- 
monic 2D XY model, but with T given by T(i?e). To 
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log (m-<m>/a) 



FIG. 1. Probability distribution function of the magneti- 
zation in a 2DXY model for different values of temperature: 
J/KgT = 0.75(c), 0.78(*),0.82(o),0.85(+) and 0.9(x). The 
solid line represents the universal distribution of BHP [1]. 



elsewhere. 

In conclusion, we have shown how universal scaling 
phenomena can arise in finite critical systems due to hy- 
perscaling; the observed similarity of the probability dis- 
tribution scaling function in the harmonic 2DXY and a 
closed turbulent cell is a special case, and is not generic. 
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see how the effective temperature should scale with Re, 
note that the velocity of the flow scales as RH, where R 
is the radius of the plates, and £1 is the angular velocity. 
The net mass per unit height of the cylinder scales as 
R 2 . Hence the net kinetic energy scales as R 4 £l 2 . The 
Re number is proportional to R 2 fl so that the energy 
scales as Re 2 . The number of degrees of freedom giving 
rise to this turbulent energy is proportional to Re 9 / 4 , so 
that T(Re) ~ Re- 1 / 4 . Using the fact that [i/a ~r@| 
we obtain that [ija ~ Re^ 1 ! 4 , which agrees reasonably 
with the data, although the exponent is not the optimal 
fit to all the data points, as shown in figure 2. 



A 

c 

v 



-.h 
O 



' Data 

- Best Fit: Slope -0.33 

— Predicted Slope: -0.25 



0.5 



4.5 



log 10 (Re) 



5.5 



FIG. 2. Data from J-F Pinton et al. [4] with the best fit 
and a fit to our predicted slope of -0.25. 

Finally, we mention that the dynamic universality class 
of the height fluctuations should be the 2D Edwards- 
Wilkinson model jlq ]; implications of this for the time 
dependent correlations of the power will be discussed 
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